Phase transition in the Ising model on a small-world network 
with distance-dependent interactions 
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We study the collective behavior of an Ising system on a small- world network with the interaction 
J(r) oc r~ a , where r represents the Euclidean distance between two nodes. In the case of a = 
corresponding to the uniform interaction, the system is known to possess a phase transition of the 
mean-field nature, while the system with the short-range interaction (a — > oo) does not exhibit long- 
range order at any finite temperature. Monte Carlo simulations are performed at various values of 
a, and the critical value a c beyond which the long-range order does not emerge is estimated to 
be zero. Thus concluded is the absence of a phase transition in the system with the algebraically 
decaying interaction r~ a for any nonzero positive value of a. 

PACS numbers: 89.75.Hc, 89.75.Fb, 75.10.Hk 



Small- world networks, which are intermediate between 
the local regular networks and the random networks, 
have two interesting features: high clustering, character- 
istic of regular networks, and shortpath length, typically 
observed in random networks P, Q . The characteristic 
path length I is defined to be the average of the shortest 
path lengths between any two nodes: In many complex 
networks, the behavior I ~ logiV with the network size 
N is observed, while regular lattices display I ~ O(N). 
When a dynamical system of many elements is put on a 
small-world network, information exchange between two 
elements involves only relatively short distance of the or- 
der of logiV. It has been observed in recent studies that 
the presence of shortcuts induces long-range order be- 
low the critical temperature Furthermore, nature 
of the phase transition has been unambiguously demon- 
strated to belong to the mean-field universality class. In 
those studies of the spin models on small-world networks, 
spin interactions have been assumed to be uniform, which 
is rather unrealistic in view of usual systems in nature. It 
is thus desirable to examine the general case of distance- 
dependent interactions, which has not been addressed. 

In this paper we study the Ising model on a small- 
world network, constructed from a one-dimensional lat- 
tice, with the interaction J decaying algebraically with 
distance r, i.e., J(r) cx r~ a . In such a case of the al- 
gebraically decreasing interaction, the critical behavior 
is expected to change with the exponent a: The limit 
a — > corresponds to the uniform interaction, which 
yields a mean-field transition at the finite critical temper- 
ature ||. In the opposite limit of a — * oo, the spins in- 
teract with only their nearest neighbors, resulting in the 
absence of long-range order at finite temperatures. This 
consideration suggests the existence of nontrivial criti- 
cal value a c beyond which the finite-temperature phase 
transition disappears. 

Here we construct a small- world network in the follow- 
ing way: First a one-dimensional lattice of N nodes is 
constructed with each node connected to its 2k nearest 



neighbors, where k is the local interaction range. Then 
each local edge is visited once and a random long-range 
connection (shortcut) is added with the probability P 
without removing local edges. Note the difference from 
the original Watts and Strogatz (WS) construction 0, 
where local edges are removed and reconnected to a ran- 
domly chosen node. 

The Hamiltonian for an Ising model on such a small- 
world network is given by 



H = -^2-1,, a, a,. 



(1) 



where <7j (= ±1) is the Ising spin on node i of the network. 
The distance-dependent interaction reads 



J it 



Jr.:, a if i and i are connected 
otherwise, 



(2) 



where is the geometrical distance (rather than the 
shortest path length) between the two nodes i and j on 
the underlying one-dimensional lattice. 

From the comparison of the two limits a — > and 
a — > oo, a total of three distinct regimes may be ex- 
pected 0-0: (i) mean-field-type critical behavior for a < 
a\] (ii) continuously varying critical exponents for a\ < 
a < 012', (iii) the short-range interaction regime with 
no phase transition for a > a.2- In Ref. 0, two bound- 
ary values a.\ and oti were obtained via renormalization- 
group calculations for the 0(n) model on a d-dimensional 
globally-coupled system. In particular, it is well known 
that ot\ — 3/2 and ol^ — 2 for the one-dimensional 
globally-coupled Ising model. Note here that the total 
number of connections in this one-dimensional globally- 
coupled system, given by N(N — l)/2 with N being the 
number of spins, is far larger than that in a small- world 
network, which is of the order of N. Accordingly, smaller 
boundary values of a are anticipated for the small-world 
network. 
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We have performed extensive Monte Carlo (MC) sim- 
ulations with the heat bath algorithm at various values 
of P and a. The range k = 2 was taken for convenience; 
other values of k (> 1) are not expected to give qualita- 
tively different results. Measured in the simulations are 
Binder's cumulant Q, the susceptibility, and the specific 
heat: 



U N = l- 



3[(m 2 }] 2 



C'r 



l(H 2 ) - (Hf 
T 2 N 



(3) 



(4) 



(5) 



with m = |(l/iV) J^i Here (•••} and [•••] denote 
the thermal average, taken over 5 x fO 4 MC steps af- 
ter discarding 5 x 10 4 MC steps for equilibration at each 
temperature, and the average over different network re- 
alizations, taken over 50 to fOO configurations, respec- 
tively. When P = 0, the network reduces to a simple 
one-dimensional lattice, displaying no long-range order 
at any finite temperature. For P ^ 0, on the other hand, 
the presence of long-range shortcuts deprives the system 
of the locally connected one-dimensional character. In 
particular, when the interaction is uniform (a = 0), the 
system at P ^ undergoes a finite-temperature transi- 
tion of the mean-field nature Q ■ 

To find a\ and a2, beyond which mean-field behav- 
ior and long-range order do not emerge, respectively, we 
have examined the behavior of the system with a varied 
at fixed P. At any nonzero value of a, the algebraically 
decaying interaction is cut off by the finite network size 
and consequently, the ground-state energy in general de- 
creases as N grows. For given values of k, P, and a, 
we thus normalize the coupling strength J so that the 
system has the size-independent ground-state energy per 
spin. As a is reduced, the interaction does not yet decay 
substantially at the boundary; this makes it inevitable 
to study systems of very large sizes for obtaining correct 
scaling behavior. 

Figure n presents the MC results for the network with 
P = 0.5 and a = O.f. Binder's cumulant Un, the sus- 
ceptibility x, and the specific heat Cy as functions of 
the temperature are shown in (a), (b), and (c), respec- 
tively. The absence of a unique crossing point of Binder's 
cumulant shown in Fig. ^a) and the shift of the suscep- 
tibility peak in Fig. ^b) as well as the peak position of 
the specific heat in Fig. ^c) unanimously suggest the 
absence of a phase transition at finite temperatures. At 
the temperatures where x an d C v display peaks, fluc- 
tuations are large and the correlation volume becomes 
comparable with the system size. This gives some infor- 
mation about the temperature at which critical phenom- 
ena are observed. Figure ^b) shows that as the system 
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action Jij oc r\y 



1: Ising model on the small- world network with inter- 
for P = 0.5, k = 2 and a = 0.1. (a) 
Binder's cumulant Un possesses no unique crossing point, (b) 
The peak temperature of the susceptibility \ keeps decreasing 
for large system sizes (N > 1600). (c) The specific heat C v 
does not display a unique crossing point; furthermore, peaks 
appear at lower temperatures as the system size gets large 
(N > 3200). 



size is increased, the peak temperature of \ increases 
first (N < 800), in accord with the observation in Rcf. 3, 
and begins to decrease eventually beyond N = 1600. In 
particular the peak temperature does not saturate to a 
certain value but decreases continuously, supporting the 
absence of a finite-temperature transition. The specific 
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FIG. 2: Logarithmic behavior of the susceptibility peak tem- 
perature [see Fig.^b)] with the inverse system size. Open cir- 
cles, triangles, and squares represent data points for P — 0.5 
and a = 0.1, 0.2, and 0.3, respectively, shown with the fit- 
ted lines. The larger the value of a is, the faster the peak 
temperature decreases. 



heat exhibits similar behavior in Fig.^c), although the 
system size at which the peak position starts to shift to- 
ward lower temperatures does not coincide with that for 
the susceptibility. 

To clarify the drift of the peak temperature, we show in 
Fig-Glthe peak temperatures T p of the susceptibility ver- 
sus the inverse system size 1/N, together with the fitted 
lines described by the logarithmic form: T p = a — b log N. 
The peak position thus approaches zero logarithmically 
as the system size grows. Here the values of the fitting 
parameters are given by a = 1.62 and b — 0.016 for 
a = 0.1; a = 1.67 and b = 0.042 for a = 0.2; a = 1.69 
and b = 0.064 for a = 0.3. This indicates that the peak 
temperature decreases faster as the exponent a is raised. 

As mentioned above, simulations here have been per- 
formed in the crossover region, giving rise to logarithmic 
behavior of the peak position with the system size. When 
the exponent a is small, the interaction reduces very 
slowly with the distance. Nevertheless it gets weakened 
substantially in a system of very large size. For example, 
when a = 0.1, the interaction between the two connected 
nodes of the maximum distance (i.e., half the system size) 
is only about 40% of the nearest-neighbor interaction in 
the system of size N — 51200. Consequently, long-range 
connections do not contribute significantly to the sys- 
tem to be ordered in the thermodynamic limit, and the 
absence of a finite-temperature transition is not very sur- 
prising. We have also investigated the case of a = 0.05, 
which requires simulations of even larger system size, to 
obtain similar results. These observations strongly sug- 
gest that the system does not undergo a phase transition 
at a finite temperature for any nonzero value of a, thus 
leading to a% = «2 = 0. Figure [21 shows that the conclu- 
sion cvi = cv2 = is apparently valid even for P = 1.0, 
which corresponds to the network with the largest num- 
ber of long-range connections. Accordingly, the common 
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FIG. 3: Susceptibility \ an d Binder's cumulant C/jv in the 
Ising model on a small- world network with interaction Jij (X 
r~- a for P = 1.0 and a = 0.1. (a) The peak temperature 
of the susceptibility decreases as the system size grows for 
N > 1600. (b) A unique crossing point is not observed in 
Binder's cumulant, implying the absence of a phase transition 
at a finite temperature. 



belief that small- world networks are similar to globally- 
coupled networks in the transition nature 0, 3 is valid 
only for the (ideal) uniform interaction. In the more re- 
alistic case of the interaction decaying with distance, the 
two systems display striking difference. 

In summary, we have studied the Ising model on 
small- world networks with the interaction decaying alge- 
braically with exponent a. For this, we have performed 
extensive Monte Carlo simulations on the networks with 
the probability for adding shortcuts given by P = 0.5 
and 1. In both cases, absence of a finite-temperature 
transition has been observed at any nonzero value of a. 
In particular, the absence in the latter, which is the ex- 
treme case of the small-world network with the largest 
number of long-range shortcuts, implies the validity for 
any values of P ( < 1 ) . This conclusion should also hold 
for the original WS network, where shortcuts are intro- 
duced by rewiring of the local edges, since it has less 
interactions than the network investigated in the present 
work. 
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